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In 1962, Scarf proposed cyclotron absorption as the physical mechanism explaining the high- 
frequency cutoff of nose whistlers. Using the Scarf proposal, two anisotropic electron velocity distri- 
butions for the magnetosphere are assumed and the complex refractive index is evaluated for each. 
The transverse (with respect to the ambient earth's field) velocity distribution in each case is assumed 
Max w ell i an. The first longitudinal distribution of velocities considered is Maxwellian to several mean 
square velocities and then is proportional to v~ l . The cyclotron damping term is then evaluated; 
however, upon investigation, the rate of change of whistler damping with frequency is found to be 
insufficiently rapid to agree with the observed whistler cutoff. The second velocity distribution con- 
sidered is a double-humped Maxwellian, i.e., a thermal electron distribution, and a resonant electron 
stream. This distribution allows for cyclotron absorption and, at the same time, is consistent with the 
whistler dispersion and attenuation. Arising from the analysis of the double-humped Maxwellian 
distribution is a transverse instability developed by Bell and Buneman which leads to VLF emissions. 
The relative importance of cyclotron damping vis-a-vis VLF emissions is examined and a qualitative 
explanation of some whistler-induced emissions is suggested. 

1. Introduction 

Whistlers are naturally occurring electromagnetic signals in the audio frequency range that 
result from the dispersion of lightning energy which has traveled through the outer ionosphere 
along the lines of force of the earth's magnetic field. The term nose whistler designates a par- 
ticular class of whistlers that exhibits simultaneous rising and falling tones starting at a particular 
frequency called the nose frequency. The dispersion of the energy from the lightning source results 
from the action of free electrons in the presence of a circularly polarized electromagnetic wave 
(whose sense of rotation is the same as that of the free electrons in the ambient magnetic field). 

A multipath whistler is a group of whistler traces originating from a single atmospheric that 
excites different magnetosphere paths. A typical frequency-time characteristic of a multipath 
whistler is illustrated in figure 1 of an accompanying article appearing in this journal [Guthart, 
1965]. The spectrogram shows the whistler amplitude, measured by the display intensity, as a 
function of time (abscissa) and frequency (ordinate). The vertical lines in the spectrogram indicate 
impulsive atmospherics. The horizontal lines indicate constant-frequency signals, usually power- 
line harmonics. The three strongest whistler traces of Guthart's [1965] figure 1 are identified at 
the bottom of the spectrogram by the letters A, B, and C. 

Each of the three marked traces is in reality a group of several whistler-components, a fact 
most easily seen in Trace B. The nose frequency of the leading whistler of each of the three traces 
is indicated by the frequency indicated to the left of the spectrogram. The whistler traces are 
derived from a single atmospheric source 1.78 sec prior to the arrival of the nose frequency of 
Trace A. This source is not shown in the figure, but its character is very much like many of the 
vertical lines seen in Guthart's [1965] figure 1. 

The nose-whistler upper cutoff frequency is rarely as much as 60 percent greater than the nose 
frequency. The upper cutoff frequency is most often extremely sharp (see annotated whistler of 
fig. 5b in Guthart [1965]). The damping at cutoff is of the order of 20 dB for frequency changes of 



1 This study was sponsored by the Stanford Research Institute. Some assistance was provided by the Air Force Office of Scientific Research on grant AF-AFOSR- 
783-65 to Stanford University. 

1403 

781-700 O - 65 - 2 



-n- 



= LATITUDE OF TERMINATION OF FIELD LINE 
AT EARTH'S SURFACE 




FIGURE 1. Normalized stream velocity (V / 'c) required for 
cyclotron resonance as a function of the maximum 
normalized frequency, m m . 



less than 2 percent. The damping rate is estimated by reducing the intensity of the whistler spec- 
trogram, using a precision attenuator, until it is no longer discernible from the background. In 
many instances, VLF emissions are triggered at the upper end of nose whistlers; in a few cases, 
the upper cutoff is defined by a more gradual rate of change of attenuation of whistler energy with 
frequency. It will be shown that the absorption and emission of VLF energy are manifestations 
of the same phenomenon. The more gradual rate of attenuation around cutoff of some whistlers 
can then be interpreted as the near cancellation of the decay and growth processes. 

The first proposed explanation for the relatively low upper cutoff frequency of nose whistlers 
(cold-plasma theory predicts that the upper cutoff frequency for whistler propagation is the mini- 
mum gyrofrequency along the path) was made by Smith [I960]. Smith elegantly derived an 
explanation of the nose whistler observations, based upon trapping of the whistler energy in ducts 
of enhanced ionization. Smith notes that the magnetosphere has a very low collision frequency; 
hence, the particles are constrained by the magnetic field to move along the field lines, giving 
long life to density irregularities perpendicular to the field. He then assumes that the duct prop- 
erties change very slowly in the direction of the field lines, so that he can represent the propagation 
path by a slab model in which the field is uniform and in which the electron density has a gradient 
only in a direction transverse to the field. For density gradients that are small in the space of a 
wavelength in the medium, ray theory is then applied. 

Using the ray-tracing approach, Smith found, for example, that for a density enhancement of 
5 percent, the whistler energy was trapped in the duct when the frequency was less than one-half 
the nose frequency, and for wave normals within 45° of the magnetic field. Furthermore the effec- 
tive group velocity of the trapped energy is most always within 1 percent of the value obtained as- 
suming purely longitudinal propagation (i.e., the wave normal and magnetic field are alined). 
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For a uniform duct, the upper cutoff frequency for trapped energy remains at SO percent of the gyro- 
frequency. For a duct whose properties change slowly in the direction of the field lines, the cutoff 
can exceed 50 percent of the minimum gyrofrequency and still be trapped in the duct. Smith's 
theory for nose whistler propagation is impressively supported by experimental evidence. Only 
the whistler properties at the upper cutoff frequency fail to satisfy the duct hypothesis. They fail 
to satisfy the duct hypothesis in two respects. First, the upper cutoff frequency of 50 percent 
predicted by the uniform duct theory is frequently exceeded. Second, the experimentally observed 
whistlers exhibit a sharp cutoff frequency. It appears that both these inconsistencies of the duct 
theory may possibly be made less extreme by the calculation of a complete ray tracing in nonuni- 
form ducts. Such a calculation would restrict the duct shape as well as the angular distribution 
of whistler wave normals. At the present time, there is no experimental evidence to support such 
restriction. As a result, a body of work has evolved which accepts the general results of the duct 
hypothesis, but looks for an alternative explanation of the high-frequency cutoff. 

As an alternative explanation of the high-frequency cutoff, Scarf [1962] and Tidman and Jaggi 
[1962] proposed the cyclotron absorption mechanism to account for the observed high-frequency 

cutoff of whistlers. Cyclotron absorption is a mechanism to account for the absorption of energy 
from a wave by particles in a collision-free plasma, and corresponds to Landau damping, with the 
principal difference being that the resonant particles absorb energy from the wave in a plane per- 
pendicular to the drift motion of the particles along the magnetic field. The absorption of energy 
from the wave is large when a large number of particles have streaming velocities such that the 
Doppler-shifted frequency of the wave as "seen" by the particles equals the particle gyrofrequency 
or its harmonics in magnitude and sense. In this article we shall consider only the gyrofrequency 
as the harmonics are cut off for a wave propagating along the magnetic field in the whistler mode. 
The resonant condition for the electron gyrofrequency is mathematically expressed by the condition 

a) — kv — ft = 

where co is the whistler radian frequency, k and v are respectively, the whistler wave number and 
electron velocity along the magnetic field, and H is 27T times the electron gyrofrequency. Because 
whistler frequencies always lie below the electron gyrofrequency, it is necessary for the streaming 
velocity of nonrelativistic particles to be in opposite sense to the wave velocity. 

The first analysis accounting for the upper cutoff frequency of nose whistlers in terms of 
cyclotron absorption assumed an isotropic electron velocity distribution function [Liemohn and 
Scarf, 1962]; however these distributions do not agree with the whistler observables in two respects 
[Guthart, 1963; Liemohn and Scarf, 1964]. First, these distributions do not account for the oc- 
currence of nose whistlers at high magnetic latitudes; i.e., for these distributions, nose whistlers 
traversing high magnetic latitude should not be observable. Second, the magnetosphere tempera- 
ture required by these distributions is too high. 

If cyclotron absorption is to explain the upper cutoff frequency of nose whistlers, alternate 
electron velocity distributions are required. In this article, a pair of anisotropic electron velocity 
distributions will be considered and tested for agreement with the whistler observations. Spe- 
cifically, a satisfactory magnetosphere electron velocity distribution must satisfy the following: 

(1) It must be stable. The general character of the upper cutoff frequency of nose whistlers 
is highly repeatable, hence the velocity distribution function responsible for the phenomenon must 
be ever present. 

(2) It must provide for a rapid attenuation at the whistler upper cutoff frequency. The rate 
is of the order of 20 dB for a frequency change of 2 percent. 

(3) It must provide for no measurable increase in dispersion beyond the zero-temperature 
estimate. A measurable increase in dispersion is approximately 2 percent. 
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2. Theory 

The refractive index, n p , of the whistler mode in a thermal magnetoplasma can be shown 
[Stix, 1962] to be given by 

2 i ■ gg? f -/ j 2 L_^£ ^i m 

p /neo^J-oc Jo L x (a) — kv z — II) 

where t; = electron velocity 
e = electron charge 
77i = electron mass 
co = radian frequency of the wave 
eo — free space permittivity 

k = wave number for the longitudinally propagating transverse wave 
fl =2tt (electron gyrofrequency) 



( v l \ 

fo = fo ( — ' v z \ = electron distribution function 



and the subscripts %, y, z define a rectangular coordinate system with the ambient earth's field, 
#o, alined along the Z-axis. The subscript ± refers to the perpendicular direction with respect 
to the static magnetic field. Consider the transverse velocity distribution to be Maxwellian, so 
that 



o *v g = exp — -^' 

2 J rrrn* a 2 



g(v z ) — z-directed distribution of electron velocity 



such that 



J: 



dv z g(v z ) ~ 1, 



a 2 = 2 k bT\ 772, 

T= transverse electron temperature, 

yV= electron density, 

A: fi = Boltzmann's constant. 
Substituting into (1) for/ , the refractive index becomes 

ia) 2 a 2 _ 

"* =1 -^ <e> < ,2) 

where 

0(v z ) 



<e> = ik\ dv 



\ 



-kv z — £V 

and 



a 1 o)\dv z a z / 



We shall now consider two longitudinal electron velocity distributions. 
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2.1. Cyclotron Absorption by Resonant Electrons With a v z l Distribution 

As a first trial function consider a distribution related to one recently proposed by Liemohn 
and Scarf [1964]. They have suggested a distribution for the resonant electrons (the longitudinal 
electron energy for the resonant electrons is typically between 0.2 and 2.0 keV) of gAVzY^Vz 1 so 
that [(l//Vo) (dNIdEzi^Ez* when dN/dE z is the rate of change of electron density with energy. 
The electron energy, E z , is equal to [(1/2) mvi\. The resonant electrons are those electrons 
participating in the cyclotron absorption, namely, those electrons whose drift velocity, Vo, is 
approximately given by 



V = - 



il 



For simplicity we shall assume that the longitudinal distribution, g r (v z ), has a v z ] beyond the energy 
range proposed by Liemohn and Scarf. Although this assumption is obviously not physically 
realistic, the distribution of the nonresonant electrons does not significantly affect the cyclotron 
absorption mechanism. The following calculation, which attempts to test the validity of the 
proposal for cyclotron absorption of whistler energy by resonant particles having a Vz A distribution, 
is thus unperturbed by the assumption on the distribution of nonresonant particles. For the least 
energetic electrons, we shall assume a Maxwellian distribution. Therefore, let 



' exp-^l/q) 2 , , , . . 

tjz \v z \ ^ sa where 5 is a positive integer 

it ' a 



(3) 



g(Vz) z 



exp (— s 2 ) sa 



rl/2, 



7r"a \v z \ 

\vz 



c > \vz\ ^ sa 



that 



f>* 



)dv, 



]+ 25exp-s^ = 



7T 1 



= T xfc&Mvz ~ a 2 + S 6XP i ^ (c 2 ~ s 2 a 2 ). 

J — oc 7T 



Let us further assume that s is sufficiently great so that v\ ~ a 2 , and that the apparent magneto- 
sphere electron temperature is approximately the Maxwellian value < 10 4 °K. The refractive 
index for this electron density distribution can be determined upon substitution of (3) into (2) and 
integrating. The integration of (2) is straightforward when it is remembered that the integral 



/ 



dv 



aa — v z 



has a pole at v z = — \a\a and is equal to 



I 



dv 



aa — Vz 



iri + In 



c+ aa 
aa-\- sa 



The quantity a is defined by 



a- 



o) — fl 

ka 
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and is less than zero for propagating whistlers. The refractive index squared for the distribution 
of (3) is thus given by 



(ajplo)) 2 



1 + 



X P (-5 2 ) f r 



^ [ 2(l+^Vn^-(y+ W 



2a 2 J *" 5a V* 7 2a: 2 
where y = l)/w, and [(1 — y)/a] has been substituted for {ka/oj). 
Let 



iri -\- In —5- 



(4) 



and 



When /£,,; < < n pr ~ /io, 



^/> ^;>/- ~r Ifl pj, 



,2 =1 + (gvM 2 . 

l ° y-l 



. _ a) 7r 1/2 &)yjili' exp (— i 2 ) 

' ~ 7 ""' 2«o 1 / 2 (fl-w) 1 ' 2 



(5) 



The total attenuation, L, experienced by an electromagnetic wave propagating in the whistler 
mode through the magnetosphere is given by 



/, 



L = 8.68 ki(6)ds decibels. 



(6) 



Since the total attenuation is linearly related to A;, the rate of change of attenuation with wave 
frequency is linearly related to the rate of change of A;(0°) with wave frequency. As noted previ- 
ously, the rate of change of damping at cutoff for nose whistlers is most often extremely sharp. 
The damping at cutoff is of the order of 20 dB for frequency changes about the cutoff frequency of 
2 percent. The rate of change of A; for the v~ l distribution is 



dki _ TT 1/2 copCls exp (— s 2 
d<o~ 2coj l / 2 (n-o)) 112 



2o)-il 



2co(il — 00) 



dki 

da) 



Zw 
(11 -co) 



A representative set of whistler parameters are [Liemohn and Scarf, 1962J: fli — 70.9(10 3 ) 
rad/sec, O = 61°, w = 39.2(l(F) rad/sec, k t {0°)= 1.5(10- 6 )m- 1 . For these parameters, the change 
in attenuation is less than 0.6 dB for a 1.5-percent change in wave frequency. As a consequence, 
this distribution does not satisfy a whistler observable (the rate change of attenuation about the 
cutoff frequency is too small); hence, it is rejected as a likely magnetosphere distribution. 



2.2. Cyclotron Absorption by a Double-Humped Maxwellian Distribution; Calculation o£ 

Refractive Index 

As a second trial function, consider a double-humped Maxwellian distribution, specified by 

(v-Vo) 2 



ivy- 



^r-exp (-^i) + -£- 



a 2 J ' n^ai 



a\ 



(7) 
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where 



/: 



= Ni 

X No 



g(v)dv=l+ X ~ I X<<U 



and 



/: 



v 2 g(v)dv = a 2 + x(aHVl) 



For simplicity, we shall further assume 



I 



v 2 g(v)dv ~ a 2 , 



although recent data reported by Serbu [1964] indicate that this is not always true. 

The velocity distribution of electrons has two components. The first is an ambient isotropic 
electron distribution defined by a temperature T= (ma 2 l2kn). The second is an electron stream 
of mean longitudinal velocity, F () , having a longitudinal temperature T\. 

As before, 



" 2/o> 2 



(2) 



After Stix [1962], the refractive index for the g(v) specified by (7) can be determined. As in the 
previous example, caution must be exercised when integrating (2). The square of the refractive 
index becomes 



%l=\ 



w(o)-il) w 2 ((o — kVo— CI) n,,((i i/c)a> 3 



u'j, 



0)1 exp (-of), (8) 






where a.\ = ; < — 1 . 



To obtain this solution, the integral 



a 2 J_ 



dv z 



I 



pl/2/ 



exp + 



cr 7T ll -(i\ 



t ex p 



(v z ~Vo) 2 



a\ 



(oj-kv z -£l) 



F - 



2X F 

a z a\ 



^her 



vl 



. . r dv z exp • 
F = jk_ r x ^ a 2 _ i_ 

TT ll2 ]-K co — fi — kv z a 



and 






dv z exp — 



(v-Vof 



a\ 



oj — kv z — CI 



h 7r^ 2 exp — a 2 . 



The quantities a and ot\ have been previously defined. It should be noted that (8) is also the solu- 
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tion for the refractive index squared for the following electron distribution function: 



/o 



N 



7r 3/2 a 3 



exp 



( ( - : 



- + 



X N 



a 2 7T 3/2 crj.ai 



exp 



exp 



(v-Vo) 2 



(9) 



The same degree of approximation is involved in the solution of (1) for the distribution of (9) as 
has been made previously. For the distribution of (9), the subscripts J_ and 1 refer to stream 
parameters only. As a consequence we shall hereafter consider (8) as the solution for the dis- 
tribution given by (9), and therefore consider the transverse and longitudinal stream temperatures 
as quantities independent of the ambient plasma temperature. For n }) = n pr -\-in P i and if rip >> 1, 
then, 



nl- 



o)(ii — o)) 



1 + 






tt/;0, 



(10a) 
(10b) 



XTT 



1/2 



(co- 



-/cFo-a)^+a 
1 1 



(O-c 



2(a l /c)a) 2 



exp 



gy-kVo-a 
O)n p o{a x lc) 



(11) 



Equation (10) specifies the real part of the refractive index as a function of the plasma and wave 
parameters. It is seen that the effect of the stream is to perturb the isotropic ambient plasma 
solution. 

Guthart [1965] has examined the dispersion of 23 nose whistlers and found the whistler dis- 
persion to be satisfied by the zero-temperature dispersion formulae; i.e., the measured whistler 
dispersion exhibits no temperature effects. (The accuracy of the dispersion measurement, which 
determines the degree to which temperature effects can be observed, is 2 percent.) If the double- 
humped Maxwellian distribution is an allowed one, then the proposed stream parameters should 
generate no observable increase in dispersion. 

Equation (11) specifies the imaginary part of the refractive index as a function of the wave 
and plasma parameters. It is seen that the imaginary part of the refractive index is primarily a 
function of the wave and stream parameters; i.e., this term becomes significant only when the 
stream velocity becomes equal to the Doppler-shifted wave velocity. The expression for the 
imaginary part of the refractive index is an even more rapid function of magnetic latitude and 
frequency than the isotropic Maxwellian solution considered previously. As a consequence, each 
frequency above the cutoff frequency is absorbed at a magnetic latitude for which the Doppler- 
shifted stream velocity is resonant with the wave velocity; hence, 

Uf)^h [6(f)] fm^f^fm. 

An interesting result of the determination of n P i is gleaned upon investigation of the term 
within the brackets of (11): 



(ft>-*r -ft)^+a 
1 1 



For a propagating magnetosphere whistler, (co — kV — ft) < 0, so that if 



ft 



■>■ 



T x \(<o-kV -Cl) 
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>1, 



(12) 



then the imaginary part of the refractive index is negative, indicating wave growth representative 
of a transverse instability. The general requirement for transverse instability (that TJT\ > 1) 
was discussed by Stix [1962] and by Sudan [1963]. The physical mechanism for this growth 
was first suggested by Brice [1963], and subsequently proved by Bell and Buneman [1964] as a 
mechanism for VLF emissions. 

Bell and Buneman considered an ambient isotropic plasma permeated by a tenuous stream 
having a mean square transverse velocity spread, v\, and a delta-function longitudinal distribution 
of stream velocity, 8(v z —Vo). Bell and Buneman determined the conditions for wave growth of 
such a distribution, and then estimated the cyclotron absorption. Equation (11) more exactly 
specifies the relative importance of growth to cyclotron absorption. To obtain significant wave 
growth, (12) requires the transverse temperature, T ± , to be much greater than the longitudinal 
temperature; the exponential factor becomes sufficiently small only when (a) — kVo — fl) — ► 0. 

We shall postpone further discussion of the growing-wave solution of (11) and calculate some 
stream parameters satisfying the whistler observables as well as establishing the stability of such 
a stream. 

3. Representative Stream Parameters 

To determine an allowed set of plasma parameters, we shall assume a value for the resonant 
particle flux, and then solve (10) and (11) for the set of plasma parameters so that (1) is consistent 
with the whistler observables and the assumed particle flux, and (2) has no longitudinal instabilities. 

The flux of resonant particles participating in cyclotron damping is given by 



r vo+bv 

nv r = vg(v)dv, 

Jl„-Ar 



where A?; is a measure of the velocity spread of the resonant particles. Jackson [1960] estimates 
this velocity to be 



i.e., those particles whose relative velocities (with respect to the Doppler-shifted wave velocity) 
are insufficient to escape the potential well of the wave contribute to the damping. 

For the double-humped Maxwellian distribution, the flux about the resonant particle velocity is 



xNai f a ' + t 

J (I, 



v-v {) v 



nv r 



yNai 



r l/2 



exp 



x 2 + 



(l\ 

AzA 



-0 + 



V 



y-t- — ) exp (— y 2 )dy, 



^ 2 j sinh ?o^v + ^ xNVn j $ r^ ^ a+ ^ 



-cD 



V5(.-£ 



(13) 



trhere 



exp (—z 2 /2)dz. 



A second constraint on the stream parameters is that the velocity distribution be stable to 
longitudinal waves, i.e., the distribution should not give rise to longitudinal plasma oscillations. 
Bell and Buneman [1964] estimate the relative significance for the growth of longitudinal waves 
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vis-a-vis the Landau damping of these waves. For a stable distribution, 

— 0.69x 1/3 <l. (14) 

a i 

The foregoing stability criterion may be unduly severe. The cyclotron absorption mechanism as 
well as the Bell-Buneman transverse instability are dependent on the steepness of the distribution 
function for velocity, v being greater than the magnitude of the stream velocity, \Vq\. The lon- 
gitudinal instability depends on the parameters of the distribution function wherein v < | Fo|, i.e., 
in that part of the distribution function where dg{v)/dv > 0. As a consequence, a distribution 
can by hypothesized that has a non-Maxwellian character for v < |Fo| and a Maxwellian form for 
v > \Vo\. Such a distribution would retain all the cyclotron absorption and transverse instability 
properties embodied in (11) and, at the same time, relax the constraint of (14). 

The final constraint on the allowed distribution is imposed by the experimental observation 
that the whistler dispersion is in good agreement (within 2%) with the cold-plasma approxima- 
tion; i.e., 



Pi _ Is (ngr—n^ds 



0.02. 



Do J s rigods 

For the double-humped Maxwellian distribution in the vicinity of resonance, 

^_ d[(o(n p r — n )] xtopiv — kVoJin — (o) 1 ' 2 
(n gr -n g0 )- -^ 2a,^(fl -o) + AK ) ' (15a) 

n(j0= 2a>^(n-o)yi*' (15b) 

Assume that the resonant particle flux, nv r , is 10 7 (electrons cm~ 2 sec 1 ), the whistler field strength 
E= 10~ 3 v/m [Carpenter, private communication], and that the whistler observables are those of 
the College, Alaska, whistler [Liemohn and Scarf, 1962]: 

/V(0°) = 250 cc- 1 

fti=7.1 (10 4 )rad/sec 
0) m = 3.9 (10 4 ) rad/sec 
n p o ~ 26 
k ~ 1.5 (10- 6 ) film. 

Using (11), (14), and (15), a set of allowable stream parameters can be derived (which is not unique) 
satisfying these constraints: 

stream density ._ _ 

y = t> \ - = 10~ 

ambient density 

ai=4(10 7 ) cm/sec 

F o = 8.4(10 8 ) cm/sec. 
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A stream with these parameters moving in the specified plasma accounts for the upper cutoff of 
the College, Alaska, whistler without increasing the whistler dispersion at cutoff. 

The double-humped Maxwellian distribution of the longitudinal velocity of electrons is seen 
to be an allowed distribution from the standpoint of whistler observables and longitudinal stability. 
Each "hump" of the distribution offers a separate and distinct contribution to the observed whistler 
characteristics. The ambient, isotropic, electrons account for the observed dispersion charac- 
teristics of the whistler. The hypothesized electron stream, with its steep gradient in electron 
density for a velocity that is greater than the mean stream velocity, | J <>|. accounts for the whistler 
cutoff properties. As previously noted, the cyclotron absorption mechanism requires a steep 
density gradient for v > |^0|. but imposes no restriction on the electron density gradient in velocity 
for v < \Vo\\ this mitigates the difficulties associated with longitudinal instabilities. 

The required stream speed for cyclotron resonance at the magnetic ecjuator is 

■ ,, I _ At - Urn _ "'», l/2 (l ~ m m fl*fflc 

\Vo\— —- ^ • (16) 

where, as previously, 

m m = Wm/fti = maximum normalized cutoff frequency, 

///i = minimum gyrofrequency along the whistler path. 

c — velocity of light in free space, and 

K\ = NI/h = constant, relating electron density and gy- 
rofrequency for the gyrofrequency model 
of electron density in the magnetosphere. 

The stream will cut off all frequencies above the maximum observable whistler frequency, co,„. 
This can be seen when it is remembered that the magnetic field strength varies along a magnetic 
field line. The minimum field strength (or electron gyrofrequency) is at the magnetic equator. 
As a consequence, a stream that is in resonance with a given frequency at the magnetic equator 
will resonate with any higher frequency at some nonequatorial location. 

Equation (16) is plotted in figure 1 for K l = 2(lO i ). It is seen from (16) that if the maximum 
normalized frequency and the electron density constant K\ are invariant, the resonant stream ve- 
locity varies as the 1/2 power of the minimum gyrofrequency along the whistler path. For the 
College, Alaska, nose whistler, the stream velocity is seen to be F () = 8.4(10 9 ) cm/sec, correspond- 
ing to an electron energy of approximately 250 eV. 

Regarding the upper cutoff frequency of "Eckersley" whistlers, it should be observed that the 
upper cutoff frequency of many of these whistlers is determined by (1) the frequency content of 
the atmospheric source, and (2) the propagation path absorption of the whistler energy. It has 
been shown [Arnold and Pierce, 1964; Dennis and Pierce, 1964] that VLF atmospherics have an 
approximately inverse-square frequency dependence in the high-frequency limit. In addition, 
the absorption experienced by the whistler energy in traversing the D layer is approximately 
proportional to the square root of frequency. Both facts limit the high-frequency content of the 
"Eckersley" whistlers, particularly at low and median latitudes, so that the rate of change of 
attenuation in the vicinity of cutoff is more gradual than is usually observed in nose whistlers. 
The source spectrum and propagation path absorption are not factors in multipath nose whistler 
cutoff for the low-latitude traces demonstrate the existence of the high-frequency components. 

The variation of stream velocity with radial distance from the earth can be inferred from the 
examination of multipath nose whistlers. A representative multipath whistler group samples the 
magnetosphere through one earth radius as measured by the nose frequencies of the whistler 
traces. Accepting the hypothesis for cyclotron absorption by a double-humped Maxwellian dis- 
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tribution, we can solve (16) for stream velocity for each of the traces of the whistler group, and so 
determine the variation of stream velocity with radial distance. A preliminary analysis of a pair 
of whistler groups leads to a stream velocity, V , distribution such that 

The existence of highly energetic electron streams is well established [O'Brien et al., 1962]; 
however, recent experimental evidence indicates the existence of streaming particles having 
energies greater than 100 eV. Serbu [1964] has reported satellite measurements out to 15 earth 
radii derived from a retarding potential analyzer that measures the energy distribution of ions and 
electrons up to 100 eV. The instrument operates in such a manner as to measure the flux into 
the analyzer of ions or electrons whose energies are greater than a given retarding potential. 
Serbu's measurements indicate the existence of electron streams having energies greater than 
100 eV. His measurements further indicate the absence of particles in the 20 to 100 eV range. 
Serbu further notes that the electrons exhibited thermal energies for geocentric distances less than 
4.5 earth radii, and that the average electron energy then increases gradually to values above 100 
eV at 8 earth radii. This would suggest that our assumption of /_* v 2 g(v)dv ~ a 2 is not appropri- 
ate. In any case, the fact does not modify the proposed analysis, which can be adapted to account 
for any observable energy. 

In addition to Serbu's measurements, Sharp, et al., [1964] have measured the energy influx 
of particles (at 300 km) precipitating on the auroral zones by particle detectors aboard orbiting 
satellites. The integrated energies of precipitating electrons of energy greater than 80 eV were 
sampled. Electron streams as soft as a few hundred eV were observed. 

Throughout this analysis we have assumed that the whistler wave normal is parallel to the 
magnetic field. When the wave normal is not alined with the field, the cyclotron resonance 
with a stream of velocity V occurs when 

| a)-k z V -£l\^0, 

where k z is the component of the wave number along the magnetic field. As a consequence, non- 
ducted whistlers would experience cyclotron absorption. In fact, Bell [1964] has recently extended 
his proposal for VLF emissions to account for nonzero wave-normal angles so that the cyclotron 
absorption and emission mechanisms are only quantitatively complicated by the introduction of 
nonzero wave normals. 

4. Whistler-Triggered VLF Emissions 

Many of the nose whistler sonograms studied to date give evidence of an emission triggered 
at the end of the whistler as, for example, in Guthart's [1965] figure 1. We can obtain a qualita- 
tive interpretation of this phenomenon by considering the previous results. It has been noted 
that 

(1) The mechanism of cyclotron absorption couples energy to the orbiting electrons. This 
energy is coupled into the plane perpendicular to the electron's drift motion along the magnetic 
field; i.e., the transverse electrjn velocity is increased. 

(2) The transverse electron stream temperature, J ± , must exceed the longitudingal stream 
temperature, T u by the factor (ft/ \co — kVo — ft |) for the Bell-Buneman transverse instability 
mechanism to dominate the cyclotron damping. 

Consider a magnetosphere having the electron stream superimposed on the ambient, thermal 
plasma. In the undisturbed state, the stream parameters are such that 

T 

1 (a)-kVo-n) 
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A whistler traversing this medium suffers cyclotron damping, surrendering energy to the normal 
component of electron motion until 



thereby triggering an emission. The electrons now return the energy previously absorbed by 
cyclotron damping, leaving the medium in the steady-state condition that existed prior to the 
whistler traversal. 

5. Conclusions 

An anisotropic particle distribution is suggested to explain the high-frequency cutoff of nose 
whistlers. The distribution consists of an ambient electron stream drifting along the magnetic 
field superimposed on an isotropic, thermal plasma. This distribution is consistent with the 
whistler observables while at the same time accounting for some VLF emissions. The suggested 
distribution is further consistent with recent magnetosphere flux measurements. 



The author deeply appreciates the guidance, insight, and stimulating counsel of R. A. Helliwell 
of Stanford University, during the course of this work, and the interest of O. Buneman of Stanford 
University. 
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